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I study the effects of particle-hole imbalance on the exciton superfluid formed in a topological 
insulator thin-film and obtain the mean-field phase diagram. At finite imbalance a spatially modu- 
lated condensate is formed, akin to the Fulde-Ferrel-Larkin-Ovchinikov state in a superconductor, 
which preempts a first-order transition from the uniform condensate to the normal state at low 
temperatures. The imbalance can be tuned by changing the chemical potential at the two surfaces 
separately or, alternatively, by an asymmetric application of Zeeman fields at constant chemical 
potential. A vortex in the condensate carries a precisely fractional charge half of that of an electron. 
Possible experimental signatures for realistic parameters are discussed. 



Introduction. — The Cooper pairs in superconductors 
known so far have zero center-of-mass momentum and, 
consequently, the order parameter is spatially uniform. 
In a singlet superconductor the two electrons in a pair 
carry opposite spin and superconductivity may be un- 
derstood as an instability of their nested Fermi surfaces. 
The nesting is lost when the Fermi surfaces are different 
in size, e.g. due to the Zeeman interaction in a magnetic 
field. But the instability may still be recovered by trans- 
lating the Fermi surfaces in momentum space so that 
they are partially nested. The pairs will then carry a fi- 
nite center-of-mass momentum equal to the momentum 
shift and the order parameter will be spatial modulated, 
as first suggested theoretically by Fulde and Ferrel [T] 
and, separately, by Larkin and Ovchinikov [2] (FFLO). 

No conclusive experimental evidence for the FFLO 
state has been found so far f3] . Recently, there has been a 
new effort [4^ in realizing the FFLO state in imbalanced 
two-species Fermi gases with magnetically-tuned inter- 
action between the two species. A different approach 
is via electron-hole bilayers where a neutral exciton su- 
perfluid is formed by the Coulomb interaction between 
the two layers. A double-layer graphene structure where 
electrons and holes are hosted on opposite surfaces in an 
external electric field is a candidate with the potential 
of a high critical temperature [5J |^. Also, a thin film 
of a strong topological insulator (STI) was argued [Tj to 
realize a novel form of the exciton superfluid, dubbed 
topological exciton condensate (TEC), where the special 
topology of the bulk results in fractionally charged vor- 
tices and protection against weak disorder [5]. Recent 
material improvements ^9 bode well for the experimen- 
tal realization of TEC. On the theoretical side, the real- 
ization of monopoles by vortices [TO] , their effective the- 
ory [TT , and the effects of screening in orbitally coupled 
magnetic field [121 have been studied. But these studies 
only consider a uniform condensate, requiring the mean 
chemical potential to be fine-tuned to the Dirac node 
(hereafter set as zero of energy) to obtain particle-hole 
balance. 

In this paper, I study the effect of electron-hole imbal- 
ance on the TEC. I show that a FFLO state is realized 




FIG. I. The phase diagram. The axes are the interaction 
strength g, temperature T, and the mean chemical potential 
of the two surfaces fi. The solid (dashed) lines indicate a 
continuous (first-order) phase transition. The spatially mod- 
ulated state is labeled FFLO and shaded blue. The dotted 
line r* (/i) separates the continuous and first-order transitions 
at finite temperature, shown by the thick green line. The in- 
sets show projections on constant T (top) and g (bottom) 
planes. 



for finite electron-hole imbalance, which can be achieved 
simply by tuning the chemical potential away from zero, 
say, by a back gate. The gap equation for the FFLO 
state in a STI thin-film is somewhat different from the 
superconducting case or the semiconductor bilayers due 
to the spinor structure of the Dirac dispersion. How- 
ever, for small imbalance, a similar phenomenology is 
obtained. I find a continuous transition from the uni- 
form to the FFLO condensate with the wave vector of 
spatial modulations set by the particle-hole imbalance 
at low temperatures, preempting a first-order transition 
to the normal state. The results are summarized in the 
phase diagram shown in Fig.[T] I also show that a vortex 
in the FFLO condensate also carries a fractional charge. 
This study suggests that the FFLO condensate is gener- 
ically present when the electrons and holes on the two 
surfaces are slightly imbalanced. Therefore, the topolog- 
ical insulator thin-film offers an exciting playground with 



several controllable knobs (separate chemical potentials 
at the surfaces, film thickness, and Zeeman fields) to ex- 
plore excitonic condensates and fractionalization in the 
absence of complications arising from the orbital effects 
of the magnetic field. 

Generalized mean-field gap equation. — In order to de- 
scribe an inhomogeneous condensate, I define a 4-spinor 

Vl'^(k)= (^^l(k+ip),^2(k-^p)), 

where Tpa is the 2-spinor at surface a — 1,2 and p 
is the wave vector of the spatial modulations of the 
condensate order parameter. Then, the noninteracting 
Hamiltonian of the two surfaces of the STI is Hq — 
Ek*p(k)^/io*p(k), with 

ho = vpTzCr ■ k + ^vpcr ■ p - AfiT^ - p., (1) 

where t and <t are Pauli matrices acting at the surface 
and spin space, respectively, Afj, = — ^2) and fi = 
-I- fi2) with the chemical potential fj,a at surface a. 
I have set h = 1. In a continuum description the sums 
are replaced by integrals appropriately. 

The Coulomb interaction can be written as a sum of 
an intralayer and an interlayer component. I will not 
consider the intralyer interaction explicitly, assuming its 
effects are taken into account by renormalizing the Fermi 
velocity. The interlayer Coulomb interaction is 

^= ^E5q'^i(q)^2(-q), (2) 
q 

where — ipaipa is the electron density in surface a, 
N is the number of sites, and = 2iTe^ eyipi^^qd) / eq by 
unit cell area, with d the distance between the layers and 
£ the dielectric constant of the intermediate medium. 

I will now derive the mean-field Hamiltonian describing 
the exciton condensate. The order parameter is 

TOk,p = ^^.gq-k(V^i(q + p)4(q))- (3) 
q 

The mean-field Hamiltonian is then i?MF = 
Eo + Ek*^(k)/i(k)*p(k) where 

Tr [mp(r)(7~^(r)m^(r)] is the kinetic energy as- 
sociated with the order parameter (the trace is over 
spin) with mp(r) oc (7(r) ^^(i/)i(x)i/'J(x — r))e~*P'^. The 
reduced Hamiltonian h is 

/i = /lo + R.e(mk,p)T^ - Im(mk,p)rj,, (4) 

The gap equation is found by minimizing the free en- 
ergy F = -Tlog(e"^*^^/'^) with respect to mk.p. As- 
suming mk,p oc 1 [7j in spin space, the result at T = 
is 

1 >p' (-gq^,p-eq,P + 5^fP^)gq-k"^q.p 
™k.P - 4^2^ £;q.,p(i?2^.p - e2_p) + „^A/ip • q ■ 

(5) 
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FIG. 2. The spectrum of energy (vertical) vs. momentum 
(horizontal). The black (gray) solid lines show the spectrum 
for finite (zero) m. The dashed (green) line shows the chem- 
ical potential at the Dirac nodes for reference, (a) p = 0; 
the particle and hole Fermi surfaces are fully nested, (b) 
p 7^ 0. The two Dirac cones in each surface are shifted by 
±|p; the solid (green) horizontal line marks the chemical po- 
tential p. = |i;f|p| for which partial nesting between particle 
and hole Fermi surfaces is obtained, (c) B, AS 7^ and p = 0; 
the solid (green) horizontal line marks the chemical potential 
Jl — BAB/Afi where particle-hole symmetry is obtained. 

Here -Eqs,p is the eigenenergy of the reduced Hamilto- 
nian Q with momentum q and spin/surface index s, 
^q.p = '^F(q^ + iP^) + + l"tq,pP, and the sum is 
restricted to occupied states with -Eqs.p < fl. The ana- 
lytical form of Eqs.p^o is complicated, but it is straight- 
forward to visualize, as shown in Fig. [2]Ja) and (b). 

In the following, I will assume for simplicity that 
5q = g = const, for |q| < A and some cutoff A. Then 
™k,p = nip independent of k. Passing to finite tempera- 
ture by introducing the Fcrmi-Dirac distribution function 
npiz) = (1 + e^^'^)^^ , the gap equation reads 

Ec^s,p{Els.p-^ls^p)+VFAlip-Cl 

(6) 

This is one of the main results of this paper. 

First, let us note that for p we recover the gap 
equation in Ref. [?]• Second, I note that the functional 
form of the gap equation is different when p ^ 0. This 
is caused by the distinct spinor structure of the surface 
Hamiltonian and is in contrast to the superconducting 
case where the FFLO state is governed by the same 
gap equation as the uniform state. However, for small 
wfIpI/too (where mg is the uniform condensate at T = 0) 
the deviations are small and we obtain a similar phe- 
nomenology. 

The change in the free energy by the condensate may 
be written as [13] 

f" da 

AF = F(m) - F(0) = - / — ^dm, (7) 
where g = g{m, fl, T) is calculated from the gap equation. 
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FIG. 3. Uniform solution (p = and units are set hy vp = 1 
and jqj < n) for Afi = 0.5 and A = 1. (a) The surface shows 
the solution to the gap equation (|6| at T = 0. Clearly the 
minimum of g occurs at m — fl. (b) The dependence of g 
on T for fi — 0.01. The values of temperature shown are 
{To = 0.0005, 0.001, 0.002, 0.003, 0.004, T* = 0.005}. 
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Therefore, a positive (negative) slope dg/dm lowers (in- 
creases) the free energy. 

Uniform condensate. — The dependence of g{m,pL,T) 
on different parameters for the uniform condensate (p = 
0) is shown in Fig. [s] At T = the minimum, (jmin, 
clearly happens at to = /I. At constant /2, the free energy 
Ai^ in Eq. ([t]) is lowered (increased) by g{m) for m > ft 
{m < jl). As a result, for a given 5, there is a finite value 
ft = p.* where these contributions cancel each other and 
there is a first-order transition from m = to* > p.* to a 
normal state. 

At r > the dip in g{m,'fi,T) is rounded and ^min 
increases with T. The upturn for to < /2 is also rounded 
and the value g{0,fl,T) initially decreases with increas- 
ing T. As a result, gnun{p,T) — ?> g{0,fl,T) continuously 
as T ^ T*{fl). For T > T*{fl), g increases mono- 
tonically. For a given g = go, let us denote by Tdfi) 
the temperature at which g(0,/i,T) goes from below to 
above go. If T*{fi) < Tdp), the transition to the nor- 
mal state will be continuous. However, if T*(/i) > Tc{fl), 
the transition will be first-order, since before reaching Tc 
the gain in free energy will be depleted by the upturn 
in g{m) at some temperature Tci{p) < T*{fi). We can 
also see this as follows: there is a value fi*{go)) where 
5* (A) = 9{^,P;T*{fl)) = gmin{p,T*{fl)) = go- For 
/2 < /i* , we have g* ( ft) < go and therefore the tran- 
sition will be continuous at some Tc{fl) > T*{jl). For 
p, > fl* , however, g*{p) > go and the transition will turn 
first order at Tdip) < T*{p). See Fig.[l] 

FFLO condensate. — The typical g(m-p,p,T 0) for 
fixed p and various values of p is plotted in Fig. |4][a). 
The corresponding free energy change /S.F is shown in 
Fig. |4|^b), which shows that the first-order transition 
from the uniform condensate to the normal state is pre- 
empted by a continuous transition to the FFLO state 
with iui?|p*| = p. This condition ensures the largest 
partial nesting between the electron and hole Fermi sur- 
faces, as shown in Fig. [2]Jb). 

At finite T, the wave vector p* associated with the 
spatial modulations of the FFLO condensate acquires 



FIG. 4. FFLO solution for p = 0.015 (A/x, A and the units are 
as in Fig.jsJ (a) g{m, p, T) and (b) The change in free energy 
AF for finite but small T = 0.00005 compared to the gap and 
four values of p = (pa.,0), = {0.001, 0.01, p* = 0.03,0.05} 
(1 through 4) with p* — 2p marked by an arrow, (c) and 
(d) show the evolution of g and AF with temperature T = 
{0.001, 0.003, 0.005, T* = 0.007} (1 through 4). Note that 
p*(r) with T -s- T*. 

temperature dependence. As seen in Fig. |4]^c) and (d), 
|p*(T)| decreases and eventually vanishes as T — > T*. 

Overall, we obtain the phase diagram shown in Fig. [T] 
where the FFLO state is shaded. I note that this anal- 
ysis is performed assuming an isotropic dispersion. As a 
result the states with the same |p| are degenerate. The 
ground state is then a superposition of such states, cho- 
sen by the residual interactions in the system or the 
anisotropy of the dispersion at higher momenta, and 
could be quite complicated. 

It was previously shown f7| that a vortex in the uniform 
condensate binds a fractional charge due, at the mean- 
field level, to a zero-energy bound state, {h + p)'i/o — 0, 
at the vortex core f6] . The spectrum is symmetric for all 
p since T{h + p)T^ = —{h -\- p), where F = TyO^yK and 
K is the complex conjugation. (Note that in real space 
k — )■ — iV in Eq. ([T]) but p is a fixed parameter.) That is, 
all states at nonzero energies must come in pairs. So, the 
single zero-energy state found for p = persists when 
p 7^ and a vortex in the FFLO condensate also binds a 
fractional charge ^e. 

Effects of Zeeman field. — A Zeeman field normal to 
surface a enters the Hamiltonian via the term Bacr^, 
where is the Zeeman energy at the surface a. This 
can be absorbed into Ec^. ([T]) by replacing k and p with 
the 3- vectors k — (k, AB/vp) and p = {p, 2B/vp) where 
AB = ^{Bi- B2) and B = \{Bi + B2). With this re- 
placement, the gap equation in Eq. ([6]) remains the same. 

In the uniform case, due to the linear density of states 
of the Dirac cone, particle-hole balance is restored for 

BAB = pA^i. (8) 
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See Fig. ^c) for a visual representation. Therefore, in 
principle, the transition to the uniform and/or the FFLO 
condensates can also be fine-tuned by an asymmetric 
Zeeman field between the surfaces. For instance with 
B = AB = 10~^ meV, and ~ 1-10 meV one can 
achieve jl/A^ - IQ-'^-lO-'^. 

Discussion. — Let us make some estimates of the pa- 
rameters. In the prototypical topological insulator 
Bi2Se3 the bulk insulating behavior is obtained when the 
thin film is a few quintuple layers, say, d ~ 10 nm. As- 
suming e ~ 30 we find g ^ 5 meV. The cutoff A ~ 1 /d, so 
taking vp ^ 1 eVA, we have ufA ~ 10 meV. Using A/x ^ 
5 meV we find 7 mo = ^/uMk^e-^^^'/s^^ - 0.1 meV. 
We should expect the FFLO state to appear when fl 
is near its value for the metastable first-order transi- 
tion. We can estimate this from the critical interaction 
strength |7] as /2 < A^e~^^^/^'^ ~ 0.1 meV, correspond- 
ing to a wavelength of modulations Afflo ^ 500 nm in 
the visible range. 

I will now turn to possible experimental signatures. 
Since the vortices are charged, the vortex flow is ac- 
companied by a charge current, which could be detected 
by various transport or noise measurements. Vortices 
would flow in response to a temperature gradient, and 
would therefore contribute to the thermopower. A vor- 
tex also has a spin texture which couples to asymmetric 
Zeeman field and can therefore be pinned by it. A vortex- 
antivortex pair has two zero-energy states and behaves as 
a two level system. Moreover, a pair carries an integral 
charge and has no vorticity. Therefore, the interaction 
between the pairs is Coulomb, and they must be decon- 
fined at T > 0, contributing to heat transport. 

The exciton condensate can sustain counter propagat- 
ing superflow in the two surfaces. Associated with this 
superflow is a large Coulomb drag |14j . whereby electric 
field or current flowing in one surface induces a counter- 
flow current in the other. The spatial profile of the cur- 
rent can distinguish the uniform and FFLO condensates. 
Recent improvements 15J in spatial resolution of local 
probes of current may be used to detect and characterize 
the structure of the superposition the FFLO condensate. 
Such high-resolution local probes could also be useful in 
detecting the counterflow in the two surfaces. 

The spatial modulations would also have optical sig- 
natures. For instance, for an FFLO condensate with a 
sinusoidal modulation (similar to the Larkin-Ovchinikov 
state in a superconductor), optical scattering experi- 
ments should be able to see changes in the reflectivity 
or refraction index of the surface [W as the wavelength 
of the incoming light is tuned through the wavelength of 
the sinusoidal modulations. 

I note that the mean-field study presented here 
tends to overestimate the critical temperature in a two- 
dimensional system, which is instead determined by a 
Kosterlitz-Thouless transition. But it is expected to cap- 



ture the qualitative structure of the phase diagram. Go- 
ing beyond the mean-field approximation, determining 
the structure of the FFLO superposition, and the effects 
of disorder are interesting problems for future studies. 

I acknowledge useful discussions with R. Budakian, H. 
Fertig, L. Li and M. Parish, and the hospitality of the 
Aspen Center for Physics where parts of this work were 
performed. This research was supported by the College 
of Arts and Sciences at Indiana University, Bloomington. 

Note added. — After this work was completed, I came 
across Ref. 171 that suggests a FFLO state may be real- 
ized in graphene double layer structures. However, the 
gap equation, the extended phase diagram, and the ex- 
perimental signatures are not discussed by these authors. 
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